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Abstract: 

This  paper  deals  with  fluid  flow  dynamics  which  may  be  Hamiltonian  in  nature  and  yet  chaotic.  Here  we 
deal  with  sympletic  invariance,  canonical  transformations  and  stability  of  such  Hamiltonian  flows.  As  a 
collection  of  points  move  along,  it  carries  along  and  distorts  its  own  neighbourhood.  This  in  turn  affects 
the  stability  of  such  flows. 


Introduction: 

There  may  be  many  settings  of  physical  interest,  where  dynamics  is  reversible,  such  as  finite- 
dimensional  Hamiltonian  flows.  In  such  transactions,  the  family  of  evolution  maps  fl  form  a  group.  The 
evolution  rule  ft  is  a  family  of  mappings  of  strips  of  transport  B,  that  we  may  consider,  such  that, 

1)  /°(x)  =  x 

2)  /'[/*'(*)]  =/t+t'(*) 

3)  (x,  t)  -*  ft{x')  from  B  x  R  into  B  is  continuous; 
where  t  represents  a  time  interval  and  t  £  R. 

For  infinitesimal  times,  we  may  write  the  trajectory  of  a  given  transaction  as, 

X(f  +  T)  =  /t+T(x0) 

=  f[f(x0,t),r] 


The  time  derivative  of  this  trajectory  at  point  x(£)  is, 

dx 
dx 


=  dTf[f(x0,t),r]\T=0  =  x(t) 

T  =  0 


The  vector  field  is  a  generalized  velocity  field, 

x(t)  =  v(x)  ...(1) 

If  xq  represents  an  equilibrium  point,  the  trajectory  remains  stuck  at  xq  forever.  Otherwise,  the 

trajectory  passing  through  x0  at  time  t  =  0  may  be  obtained  by, 

t 

x(t)  =  /f(x0)  =  x0  +  I  dx  v[x(t)],x(0)  =  x0  ...  (2) 


The  Euler  integrator,  which  advances  the  trajectory  by  St  xvelocity  at  each  time  step  is, 

Xi  =Xi  +  Vi(x)Sr. 

This  may  be  used  to  integrate  the  equations  of  the  dynamics  concerned. 


Hamiltonian  Chaotic  Dynamics: 

For  a  Hamiltonian  H(q,  p)and  equations  of  motion 

dH  dH 

dpi  dqt 

With  D  degrees  of  freedom, 

x  =  (q,p), 

(J  =  0h,<?2.<73. -.^d). 

P  =  (Pi,P2.P3.-.Pd)- 

The  value  of  the  Hamiltonian  function  at  the  state  space  point  x  =  (q,  p)  is  constant  along  the  trajectory 

x(t).  Thus  the  energy  along  the  trajectory  x(t)  is  constant, 

d  dH  dH  dH  dH      dH  dH 

dt  dqt  dpi  dqtdpi      dptdqt 

The  trajectories  therefore  lie  on  surfaces  of  constant  energy  or  level  sets  of  the  Hamiltonian 
[(q,p):H(q,p)=E]. 

Given  a  smooth  function  g(x),  the  standard  map  is, 

xn+l  =  xn   '   yn+1 

yn+i  =  yn  +  g(.xn). 

This  is  an  area-preserving  map.  The  corresponding  nth  iterate  Jacobian  matrix  is, 

K=n 

Let  M  =  1  as  the  map  preserves  areas  and  also  B  is  symplectic. 

The  standard  map  corresponds  to  the  choice  g(x)  =  k/(2n  sin(27rx)  ).  When  k  =  0,yn+1  =  yn  =  y0, 
so  that  angular  momentum  is  conserved,  and  the  angle  x  rotates  with  uniform  velocity, 

*n+i  =  xn  +  y0  =  x0  +  (n+  l)y0  ...  (3) 

The  standard  map  provides  a  stroboscopic  view  of  the  flow  generated  by  a  time  dependent  Hamiltonian, 


H(x,y,t)=-y2  +  G(x)S1(t); 
where  S1  denotes  the  periodic  delta  function, 

oo 

St(t)  =    V    S(t-m) 

m=-oo 

and        G'(x)  =  —  g(x)         —  (4) 

A  complete  description  of  the  dynamics  for  arbitrary  values  of  the  nonlinear  parameter  k  is  fairly 
complex.  When  K  is  sufficiently  large,  single  trajectories  wander  erratically  on  a  large  fraction  of  the 
phase  space. 


Stability  In  Hamiltonian  Flows: 

The  equations  of  motion  for  a  time-independent  D-degrees  of  freedom  Hamiltonian  can  be  written  as, 

/  o      —  A  d 

it  =  (OijHjix),  (»  =  {_!      Q),  Hj(x)  =  —H(x); 


Xi 


where  x  =  (q,  p)  £  B  is  a  phase  space  point.  HK  =  dKH  is  the  column  vector  of  partial  derivatives  of  H, 
I  is  the  [D  x  D]  unit  matrix  and  a)  the  [2D  x  2D]  symplectic  form. 

o)T  =  — O),  a)2  =  —  1 

The  evolution  of  ]l  is  determined  by  the  stability  matrix  A, 

-/(x)  =  AWfix), 

Aij(x)  =  G)ikHkj(x); 

where  the  matrix  of  second  derivations  Hkn  =  dkdnH  is  the  Hessian  matrix.  For  symmetry  of  Hkn  , 

ATa)  +  (oA  =  0 

This  is  the  defining  property  for  infinitesimal  generators  of  symplectic  or  canonical  transformations, 
which  leave  the  symplectic  form  a)  invariant. 

For  an  equilibrium  point  xq  the  stability  matrix  A  is  constant.  The  characteristic  polynomial  of  A  for  an 
equilibrium  xq  satisfies, 

det(A  -  XI)  =  det[a>-1G4  -  AI)co 


=  det[— a)Aa)  —  AI] 

=  det(AT  +  AI)  =  det(A  +  AI) 


The  symplectic  invariance  implies  that  if  A  is  an  eigenvalues,  then  -A,  A*  and  —A*  are  also  eigenvalues. 

A  Floquet  multiplier  A  =  A(x0,  t)  associated  to  a  trajectory  is  an  eigenvalues  of  the  Jacobian  matrix/ 
and  it  satisfies 

det(/  -  A/)  =  det(JT  -  A/)  =  det(-a)JTa)  -  A/) 

=  detCT1)  det(/  -  A/) 

=  A2Ddet(J-A~1I)        ...(5) 

This  is  because,/-1  =  —o)JTa),J  being  symplectic.  If  A  is  an  eigenvalue  of/  so  are  -,  A*  and  — .  Real 

eigenvalues  always  come  paired  as  A,-.  The  complex  eigenvalues  come  in  pairs  A,  A*,  |A|  =  1,  or  in 

1  i 

loxodromic  quartets  A,-,  A*  and  — . 

M  A  A* 

For  a  trajectory  originating  near  x0  =  x(0)  with  an  initial  infinitesimal  displacement  <5x(0),  the 
flow  transports  the  displacement  <5x(t)  along  the  trajectory  x(x0,  t)  =  ft{xQ). 

This  infinitesimal  displacement  is  transported  along  the  trajectory  x(x0,  t),  with  time  variation  given 

by, 


l^(*°'t)  =  ZSW 


;       -] 


Sxj(x0,t)  ...(6) 

X=X\Xn,Z) 


The  system  of  linear  equations  of  variations  for  the  displacement  of  the  infinitesimally  close 
neighbour  x  +  Sx  follows  from  the  flow  equations  by  Taylor  expansion  to  linear  order, 


xt  +  6xt  =  vt{x  +  Sx)  *  vt(x)  +  2_.^5xJ  -  (7) 

j 


/—i  dXj 


Taking  these  together,  the  set  of  equations  governing  the  dynamics  in  the  tangent  bundle 
(x,  Sx)e  TB  obtained  by  adjoining  the  d-dimensional  tangent  space  Sx  e  TBX  to  every  point 
x  e  B  in  the  d-dimensional  state  space  B  c  RD,  may  be  written  as, 

xt  =  ^(x),5xt  =  2^Aij  (x)5xj 
j 

and  the  stability  matrix  or  the  velocity  gradients  matrix  may  be  written  as, 

A"M=^r    -(8) 

Equation  (8)  describes  the  instantaneous  rate  of  shearing  of  the  infinitesimal  neighbourhood  of  x(t)  by 
the  flow. 

Description  of  neighbourhoods  of  equilibria  and  periodic  orbits  is  afforded  by  projection  operators, 


nM  -  A^I 

where  matrix  Mis  either  equilibrium  stability  matrix  A,  or  periodic  Jacobian  matrix/ restricted  to  a 
Poincare  section.  For  each  distinct  eigenvalue  Aw  of  M,  the  columns/rows  of  Pt  are, 

(M  -  A^)Pj  =  Pj{M  -  A^l)  =  0 

are  the  right/left  eigenvectors  e^K\  e^  of  M  which  span  the  corresponding  linearized  subspace, 
provided  M  is  not  of  Jordan  type.  For  determining  the  eigenvalues  of  the  finite  time  local  deformation  J*- 
for  a  general  nonlinear  smooth  flow,  the  Jacobian  matrix  may  be  computed  by  integrating  the  equations 
of  variations, 

x(t)  =/t(x0),5x(x0,t)  =;t(x0)5x(x0,0) 


For  equilibrium  point  xq, 


f(xq)  =  eA«t,Aq=A(xq)  ...(10) 


Conclusion: 

In  Hamiltonian  dynamics  of  real  physical  entities,  the  difference  between  support  towards  a  flow  and 
resistance  against  it  not  only  causes  the  transactions  to  take  place,  but  also  distorts  the  neighbourhood 
of  the  flow.  For  sufficiently  large  values  of  K,  a  given  Hamiltonian  flow  may  turn  into  a  highly  chaotic 
one.  In  accordance  with  the  Laws  of  Physical  Transactions  as  per  the  Theory  of  Abstraction,  such  flows 
will  generate  disturbances  in  all  parts  of  the  vicinity  of  the  transactions.  These  disturbances  will  tend  to 
smooth  out  equally  in  every  possible  direction,  causing  secondary  disturbances  against  stability. 


